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SUMMARY 


This paper encompasses that part of error analysis which deals with errors 
resulting from the instrumentation used in measuring pertinent wind-tunnel 
parameters. The parameters, for this discussion and analysis, are limited to 
those required for the wind-tunnel model measurements. The pertinent parame- 
ters, their standard deviations, and the theoretical derivation of them, are 
given. BASIC programs and plots for the standard deviations of dynamic pres- 
sure, Mach number, and Reynolds number are included for the National Transonic 
Facility. A literature search was made back to 1934 and it was found that 
information of this kind is almost nonexistent. 


INTRODUCTION 

In any successful research endeavor, it should be possible to experimen- 
tally validate the results. In addition, the same or other experimenters 
duplicating the parameters of the experiment should be able to get the same 
results within experimental error. When this cannot be achieved, it must be 
concluded that the theories or formulations that were to be validated were 
incorrect or that large errors of some kind entered into the experiment, which 
prevents duplication of results within the experimental error. 

In numerous reported wind-tunnel measurements, it is indicated or can be 
inferred that correlation of measurements between similar tests has not reached 
a level of desired consistency. The purpose of this paper is to determine the 
extent to which the accuracy and precision of the measuring instrumentation 
contributes to the total error. Though it is not dealt with herein, the two 
other types of errors are those introduced in establishing and measuring perti- 
nent tunnel parameters and those entering into the acquisition and reduction of 
data. 


Before any results can be compared, it is imperative that their accuracies 
and precisioBS be reported in a consistent manner. The major experimental 
errors are either systematic or random (accidental) . Although systematic 
errors do exist and occur, they are not statistical and can be corrected or 
eliminated. Therefore, they are not further discussed in this paper. This 
work assumes that the errors are random (accidential) , which means that they 
are variable in magnitude, follow no pattern in occurrence, and can be either 
negative or positive. General references 1 , 2, and 3 show that the arithmetic 
mean of the measured value is the one that most closely approaches the true 
value of the measured quantity. The standard deviation, which is the root- 
mean-square value of the deviations, best describes the scatter of the mea- 
sured values from the mean value. In reporting any experimentally measured 
results, therefore, the accuracy (arithmetic mean) and the standard deviation 
of the mean value a/^/n should be given or the mean and a preferably with 
n included. When grappling these results, either set of values should also be 
shown. In the following pages and in appendix A it is shown how to obtain 


the standard deviations of significant parameters of wind-tunnel measurements. 
Knowing the number of observations taken to obtain the standard deviation, the 
standard deviation of the mean value is obtained. From these and their mean 
values, it can be determined where accurate and precise measurements are more 
critical and necessary. This may require more research and development to 
advance the state of the instrument art as well as the other areas involved 
in the final wind-tunnel data accumulation. 

Since the current problem is a manifold problem, this paper deals only with 
the accuracy and precision of the measuring instrumentation being used in wind- 
tunnel testing at present, the theoretical analyses which show where the errors 
are produced, their criticality and how they combine, the development of the 
second-order strain-gauge balance interactions and how they combine, and why 
higher-order terms are ignored. Then, using the parameters for the National 
Transonic Facility (NTF) at Langley Research Center and the instrument errors 
for the instruments used in this facility in the derived formulas, the follow- 
ing standard deviations were obtained; 0(Reynolds number), 0 (Mach number), 

0 (dynamic pressure), 0 (normal- force model), 0 (axial-force model), 0(pitching- 
mcment model) , 0 (yawing-moment model), 0(rolling-mc*nent model), and subsidiary 

sigma values for items such as normal-force coefficient, side-force coefficient, 
lift and drag coefficients, etc. Programs have been written (see appendix B) 
for the Hewlett-Packard 9830A calculator for obtaining 0(Mach number), 
0(Reynolds number), and 0(dynamic pressure), and results for the NTF have been 
plotted. Assuming the values for the parameters and the instrumentation for the 
NTF remain unchanged, potential users of the NTF should find the error results 
of this report useful. If changes in the values do occur, these can readily be 
taken care of in the data section of the Hewlett-Packard program. 

Use of trade names or names of manufacturers in this report does not con- 
stitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


SYMBOLS 

A reference area 

At, local area (base) where Ap acts across base of model 

Bj^ balance output including second-order terms; for example, B-j is 

balance normal-force output; see equation (11) (i = 1 , 2, . . . , 6) 

b reference length 

axial-force coefficient 

Cab base axial-force coefficient 

Cj 3 drag coefficient 

Cdb base drag-force coefficient 
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lift coefficient 


rolling-moment coefficient 
pitching-moment coefficient 
normal-force coefficient 
yawing-moment coefficient 
pressure coefficient 
side- force coefficient 
mean geometric chord 
specific heat at constant volume 
axial- force component on model 
normal-force component on model 
side- force component on model 

jth ccanponent factor interacting to contribute to ith component; 
for example, K-j 2 is side-force component factor interacting 
to contribute first-order terms to normal force; see table I 
(i = l, 2, • • •, 6; 2, • • • f 27) 

reference length 

Mach number 

rolling moment 

pitching moment 

yawing moment 

number of observations 

static pressure 

base static pressure 

local static pressure 

stagnation pressure 

free-stream static pressure 


dynamic pressure 


free-stream dynamic pressure 
R Reynolds number 

S reference surface 

T absolute tonperature 

total tonperature 
ot angle of attack 

3 sideslip angle 

Y ratio of specific heat of tunnel gas at constant pressure to that 

at constant volume 

dynamic viscosity 

0( ) standard deviation 

There are six balance components Fjg, Fy, My^ and M 2 ) . A super- 

script ”(2)” on any of these subscripts indicates the second-order term of that 
component. For example, ^^(2) second-order term of the normal-force 

component. Any two of these components together indicate the interaction of 
the second term upon the first. For example, FjjFy is the cross term whereby 
the side- force component contributes to the normal- force output. 


ERROR ANALYSIS 

When considering the following sections, it should be understood that small 
errors in the independent variables are assumed in the resulting formulas. If 
the errors become so large that second- or higher-order terms are not negligi- 
ble, then the errors are derived by using A values and not exact differen- 
tials. It is assumed that random errors are being dealt with whose magnitudes 
are variable and whose occurrence is disordered. It is also assumed that their 
distribution is normal or Gaussian if systematic errors are removed and good 
sampling practices employed. In those equations that are applicable, the val- 
ues used for the constants are for air. Where different substances are used, 
the appropriate constants for that substance should be used. 

The general texts on errors given in the references show that the mean 
value of a measured quantity (assuming randan errors only) approaches the true 
value of that quantity. If systematic (bias) errors exist, then the mean value 
becomes displaced from its true value and results in inaccuracy. Systematic 
errors are not statistical and by definition can be eliminated when they are 
discovered. The general reference texts also show that the root-mean-square 
value of the measured quantity equals its standard deviation from the mean 
value. The mean value and the standard deviation can thus define the accuracy 
and precision of a random measurement and should always be used together when 
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this is done. If a function f(x,y,z) is dependent on the three independent 

variables x, y, and z, the standard deviation may be obtained by taking the 
exact differential of the function, e.g., 


df = — ox + — oy + — oz 
9x 3y 9z 


( 1 ) 


If df is the deviation frcxn f, then the right side of equation (1) must be 
the sum of the deviations of the three independent variables. Using superposi- 
tion and the definition of O, then 



The mean value of f is as follows: f = 


This section is a summary of the standard deviations which are valid over 
all Mach numbers and assume the use of differentials (small errors) . The deri- 
vations for these values as well as the remainder of the balance output errors 
are shown in appendix A. 

Standard deviation for Mach number (ref. 4) : 
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Standard deviation for dynamic pressure (ref. 4) : 


q = - 

2 


( 6 ) 


1/2 


a(q) = 


ryM2 


a(p) 


[ypMa (M) ]" 


(7) 


a(q) 


a(p) 


-^2 


20 (M) 
M 


n2 


1/2 


( 8 ) 


Standard deviation for Reynolds number per meter: 

lV2 


- = 4.790 X 1q4 
I 


, ,-4/7 , ,-6/7 

^Pt\ /Pt' 


Pt5/7p2/7 + 


110.3pt' 


(9) 


r 


. 1/2 


:4.790 X 104)^< 


,4/7 , ,6/7 

,S -(j) 




X (-1) 


,2/7 Tt + 220.6' 


Pt, 


Pt 


V 

) 

/ 


/rr 


- (a" 


1/2 


-f- 

7\pt 


2/7 


'p\ 


2/7 


-1 — 1 Tt + 110.3 I — 1 T<- + 110.3 


Tt^ 


Pt 5 

Tfc2 P 2 \L\Pt/ 


x4/7 , Ja/l' 

-(a . 


,-V2 


(Equation continued on next page) 
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For an interaction correction of a six-ccanponent strain-gauge balance, it 
is assumed that each component contributes second-order terms and that the 
second-order terms are small. Taking two components at a time to obtain the 

6 ! 

second-order terms (refs. 5 and 6), results in 6Co = = 15 . Fifteen 

2 ! ( 6 - 2 )! 

cross terms plus 6 first-order terms plus 6 squared terms result in 27 terms 
for each of the 6 model components obtained from the strain-gauge balances. 
These terms are listed in table I. The symbolization (see table I) of a con- 
tribution of side force, for example, on the normal force is K -|^2 times Fy 
for the first-order term and K-] 3 times Fy for the second-order term. This 
notation is more easily grasped in the follcwing matrix equations. Terms 
involving higher than second order are usually too small to be significant and 
will be ignored. The balance output column matrix is 

Normal-force output from balance 
Side-force output from balance 
Axial-force output from balance 
Pitching-moment output frcm balance 
Yawing-moment output from balance 
Rolling-moment output from balance 

^ ( 11 ) 


— j_®i^ (i — 1, 2, . . ., 6) 
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force-mcxnent equations for a six-component strain-gauge 








If we let B equal the matrix in the left c»lumn, Y equal the matrix 
in the middle c»lumn, and Z equal the matrix in the right column of equa- 
tions (12) , and let K equal the 6x6 matrix and K-| equal the 6 x 21 matrix, 
then 


B = KY + K-i Z 


Since B equals the output from the strain-gauge balances and Y equals the 
actual loads applied to the balance, solving for Y yields 

K~**KY = - K“'*KiZ 

Since K“^K, if K“^ exists, is the identity or unity matrix K”^K = D ] and 
y = - K*^K*|Z. From this, a(Y) can be obtained, which is the same as 

a(F^), a(FY), o{Fp^), a(My)r a(Mz),or a(Mx). 

Since obtaining Y is prolonged and tedious and since it, as well as 0(Y) 
is not available as yet, the development in this paper proceeds with the outputs 
of the strain-gauge balance components, matrix B and cr(B), that is, a(B-|), 
a(B2) , a(B3), a(B4), a(B5) , and aCBg) • From this it is shown how to obtain 
the extrema and thus the errors of the Y matrix. As shown in equation (13), 
a(B’|)/B'| = f (Fjq,FyfF;^,MY,Mz,Mx, the interaction coefficients, and measurement 
errors) . Since the only time F^, Fy, F;^, My, and Mx are known is 

at balance calibration, the K factors are obtained using known loads having 
known load standard deviations and using statistical methods for arriving at 
the expectation values for K and cr(K). The measured values of K are 
obtained by applying the loads singly and in combination and observing the 
effects on and the values of the outputs of interest. Once this information is 
obtained, wind-tunnel measurements can be made. For example, suppose measure- 
ments of the output of the normal component of the balance are made. If the 
loads on the model are not known, all of them may be assumed to be acting. 

Each balance component output is then multiplied by its sensitivity constant 
which was obtained during balance calibration. This yields the uncorrected 
balance outputs; i.e., CjjB-j , CyB2, ^n®5^ ^Z®6- These 

outputs are then corrected by taking the products of the interaction coeffi- 
cients and their uncorrected balance outputs and subtracting them from 
which yields a corrected B-| . Corrected B2, B3, B4, B5, and B5 are 
obtained in the same manner- The process is reiterated until the corrected 
balance outputs converge to a value which is within 1 percent of their last 
value. With these values, the standard deviations of the balance outputs can 
be found by using equation (13) in the body of the paper and equations (A21) 
to (A25) in appendix A. Once the corrected loads are determined, the calibra- 
tion curve shows the corresponding balance outputs. The positive and negative 
standard deviations can be projected onto the calibration and in turn onto the 
axis for the actual load values. For linear calibrations, these projections 
would show the load standard deviations. For nonlinear calibrations, these 
projections would show the extrema from which the errors are obtained. It is 
important to note that all of the interaction coefficients K, excluding sys- 
tematic errors, have random errors associated with their measurement, and for 
this reason they are treated as random variables because they can vary between 
the extremes set by their error limits. 
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The standard deviations of B-\ , B2, B3, B4, B5, and Bg are typified 
by equation ( 13 ). The derivation for a(B^)/Bi, in which all the K factors 
are considered to be random variables, is as follows: 


a(Bi) 


= <(1 + 2FnKi ^7 + Ki^i 3 Fy + Kij 4 Fa + Ki j gMy + Ki^igMg + Kij7%)' 


/Fn\^ 

I 

Q 

1 

\Bl j 

. J 


+ (Ki,2 + 2FyK-|^8 ^nKi ,13 + Fa^ 1,18 +%Ki ,19 


+ MzKi ^ 20 + ^ 21 ) ' 


/Fy\2 

■a(FY)' 

\®V 

_ _ 


n2 


+ (K1,3 + 2 FaKi ^9 + FijKi j4 + FyK^ j g 


+ MyKi ^ 22 + MgKi ^ 23 %Kl , 24) 


/Fa^ 

-a (Fa)- 

\Bi j 

_ Fa 


-i2 


+ (Ki^4 + 2MyKi JO 


+ FnKi J5 + FyKi j 9 + F^Ki j2 + MzKi ,25 + MxKl , 26 )‘ 


/MyN^ 

~0 (My) - 

\Bl / 

My 


+ (Ki, 5 + 2MzKi ,11 + FnKi,16 + FyKi , 20 + Fa^I , 23 + %Kl , 25 + %Ki,27)' 

i2 


/Mz\2 

■0(Mz)- 

\Bi j 

Mz 


+ (Ki,g + 2 MxKi,i2 + Fi,Ki,i 7 + FyKi , 21 + F^Ki , 24 


/Mx\^ 

-a (Mx) - 

^ /FyKl,2\2 

-a(Ki,2)- 

^ /W,3\ 

\Bl / 

% 

B1 j 

_ K1,2 

\ Bl / 


■0(Ki,3)“ 

^ /%Ki,4\^ 

-a(Ki,4)- 

^ /MzK1,5\2 

■(1(K1,5)‘ 

^ /MxK1,6V 

_ Ki,3 

V / 

_Ki,4 

Bl / 


V Bl j 


(Equation continued on next page) 


10 






X 


P<J(Ki,6)~ 

^ f%^Ki,7y 

-0(Ki j)- 

2 

■<J(K1,8)“ 


'’ \ Bi j 

_ K1,7 

\ Bl / 

_ Ki,8 


f^A^Kl , 9 ^ 

■<J(Ki, 9 )" 

\ Bt / 

_ K 1,9 




.2. 


■a(Ki jo) 

Kl,10 




■^(Kl,ll)“ 


■<^(Ki,12)' 

2 ^FnFyKij3>^2 

j3)- 

^^1,11 

J 

_ ^1,12 

\ Bi j 

_ Kl ,13 


/FNPAKnUy 

■0 (Ki,i4)~ 

2 /P#yKij5N^2 

■<J(K1,15)' 

^ /F#zKi,16\ 

\ Bi J 


\ Bi j 

_ K1,15 

\ / 


■'^(Kl ,16)' 

^ /FnMxK1,17\^ 

■Cf(Ki j7)" 

^ /FyFaK1,18\^ 

■‘^(Ki,18)' 

_K1,16 

\ Bi ) 

_ Kl j7 

^ \ Bi / 

_ Kl j8 


*i2 


/'FyMyKi jgN 
__ 


■0 (Ki j 9) 

K 1,19 


^ /FYMzKl,20^^‘ 


■0(Kl ,20) 
^1,20 


[2 /f^xKi,2iV 

n ^ T - 


-0(Ki ji)- 

2 /FaMyKi 

■0(Ki,22)‘ 

^ /"FaFyKi,23\2 

-0(Ki j3)- 

_ Kl ji 

V j 

_ Kl j2 

^ \ Bi j 

_Ki,23 


/FaMxKi, 24 V' 


■0(Ki j 4 ) 

Ki ,24 


2 /MxMzKi j 5 > 


\ 


■a(Ki J5) 


(Kl , 26 ) f /MzMxKl , 2 ?\t , 2?) 


Ki ,26 




K 1,27 


Ki ,25 
1/2 


/MyMxKij6Y 

] 


The derivations for 0(82) /Bzf 0(83) /B3, 0(84) /B4, 0(85)785, and 

0 ( 8 g)/ 8 g are given in appendix A. 


( 13 ) 
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Equations ( 12 ) are the force-moment equations for Bi , B2, B3, B4, B5, 

and Bg. At this point, since we have the parameters that define measurements 
of interest to wind-tunnel researchers, these are summarized. 

Body-axis measurements (balance axis) if balance is aligned with model 

Axial-force coefficient = 

Normal-force coefficient = Fjj/qS 

Side-force coefficient = Fy/qS 

Pitching-moment coefficient = My/qSc 

Rolling-moment coefficient = Mjj/qSb 

Yawing-moment coefficient = M2/qSb 

Angle measurements 

Sting angle 

Angle of roll 

Angle of sideslip 

Angle of attack 

Stability-axis coefficients 

Drag coefficient = f (C;^»<^»Cq3) = Cq = cos a + Cjj sin a 
Lift coefficient = f (Cjj ,01,033) “ “ *-N ^ ~ ^ 

Side-force coefficient = f(Cy,B) 

Pitching- moment coefficient = f(Cm) 

Rolling-moment coefficient = f(Ci,OL) = Cj cos ot + C^i sin a 
Yawing-moment coefficient = f(Cj^,a) = cos a - Cj sin ot 
Base axial-force coefficient = ^(Pb ” Pqo)/^S 
Base drag- force coefficient = f(C;^3,a) = Ab(Pb “ Poq) ci/qS 
Lift-drag ratio L/D 

Plot of Cl versus C3 
Cp = (p) - pj/qoo 
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Using the preceding list, the standard deviations for these coefficients 
are enumerated in the remaining equations. Where it is believed to be neces- 
sary, the full derivation for obtaining them is given in appendix A. 


o (Axial-force coefficient) 
Axial-force coefficient 



d(S) 

is self-explanatory. 

S 


a (Normal- force coefficient) 
Normal-force coefficient 


a(Side-force coefficient) 
Side-force coefficient 



( 14 ) 


( 15 ) 


( 16 ) 


a(Pitching-mc»nent coefficient) 

Pitch ing-rocxnent coefficient C„, 



+ 



( 17 ) 


a (Rolling-moment coefficient) ^(Cj) j 

"o(Mx)' 

2 

4. 

"Cf(q)“ 

2 

1 

~o(S)~ 

Rolling-mcment cx^efficient 

Mx 

T 

_ q _ 


^ s ^ 


+ 



( 18 ) 


a(Yawing-mc«nent coefficient) 

a(Cn) ( 

r 

O (M2) 

1 

r^f + 

'o(S)' 

Yawing -moment coefficient 


Mz J 

J 

_ s _ 



( 19 ) 
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The standard deviations for the angle measurements are not derived but 
measured. The technique for obtaining the angles is to use an accelerometer 
employing a null servo system. By this means a strut angle can be measured to 
a relative precision of 20 = ±0.02° and a model angle of 20 = ±0.03°. The 
acceleraneter must be isolated from vibration which has an adverse effect on 
it. Tunnel investigators would like to obtain 20 = ±0.01°. This is possible 
to obtain but is not routinely feasible at present. 




Of cxjurse, other factors contributing independently to any of these stand- 
ard deviations, such as the error in Cd and Cj, due to grit effect, internal 
friction effects, etc., must be added in. Fjj, P^, a, p^., p, and S may be 
dependent on other parameters, so their errors would be formulated on more 
fundamental parameters. This continues until the last desired quantity is 
formulated in as fundamental a system as desired. The latter is usually mass, 
length, time, temperature, and charge. This would be the suggested error- 
analysis format reconroended for use in instrument work. For wind-tunnel appli- 
cation or for any other specialized field application, the error- analysis format 
would be stated in terms of basic parameters and could be presented thusly; 


ct(Cd) = 


(sin a)a(Pu) 
qs 


”|2 r 


(cos <x)a(Fp^) 

qs 


n2 


+ (7 (a) 2 + 


Cq O (q) 


i2 


, 1/2 


(24) 


where 


'Cd o (q) 


~\2 


= <Cd 


a(p) 


2Cd a(M)l 


M 


a (Cl) 



fsin a' 


qS 


a(FA)2 + a (a) 2 + Cl^ 



(25) 


O (Cy) 

(Same as eq. (16)) (26) 

Cy 


‘^(Cm) _j 

O (My) 

2 

1 

"a(q)" 

2 

I 

"a(S)“ 

2 

4_ 

0(c) 

Cm ' 1 

« My 

T 

_ q 

T 

_ s . 

r 

_ c ^ 



(27) 
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o{Ci) 


a(Ci) 


a(Cn) = 


sin /cos a\^ /mx cos a\^ ^Mx cos 

0(a)2 + j — _]a(Mx)2+| ; ) a(q)2 + 


qSb 


qSb 


q2sb 


qSb'* 


1/2 


X a(b )2 + 


^Mx cos a' 


qbS'' 


0(S) 


(28) 


I 


O (Mj^) 

L % J 


+ 2 


0(q) 

q 


+ 2 


<7 (A) 


-i2 


+ 2 


a(b) 


i2 


+ (a tan a) 2 


0(a) 


a 


, 1/2 


0(M2> 


Mn 


+ (a cot a) 2 


0 (a) 


a 


(29) 


(Hz sin a\^ /cos a\^ . 

0 (a) 2 + ( — 0 (M2) 2 + 


qAb 


qSb I 


(Hz COS a' 
V q^Sb ) 


0(q)2 + 


/M2 COS 


qs2b 


2 2 

/M2 cos a^ ^ ^Mx cos a^ ^ /sin a\^ 


X 0(S)2 + cr(b)2 + 

\ qAb2 / 


qSb 




qSb 


/Mv sin a^ 


q2sb 


a(q)2 + 


rMx sin aN 
qs2b 


0(A) 2 + 


/Mx sin a' 

qSb2 


CT(b) 


1/2 


(30) 


CT(Cn) 


0(M2) 


-i2 


Mq 


+ a 2 (tan a + cot a) 2 


■0(a)' 

2 

+ 2 

"a(b)' 

a 


_ b _ 


vV2 


(31) 
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^ (Cab) “ 


a (Cab) 
Cab 


^ (Cdb) ~ 


Of(CDB) 

CdB 



f Pb - p„ 
^ qs 


a(Ab)2 + 


Ab(pb - pj 
q2s 


-|2 


CJ(q) 


(32) 


a(Ab) 

Ab 


“i2 


( t 

"c(Pb)' 

\Pb - Poo / 

_ ^ 


-|2 


. 2r 


(Pb - Po 


<^(Poc) 


a(q) 


-i2 


a(S) 


-.2 


)/2 


(33) 


Ab(Pb - Poo) sin a 

qs 


t Ab oos a\ /Ab c»s 

a(pb)2 + ( 

qs / \ qs / 




~\2 


Pb - Po 

V qs 


c»s a 


cJ(Ab)2 + 


Ab(Pb - Poo) cos a 
q2s 


a(q) 


'Ab (Pb - Po.^ ' 



(cos a)2a(s)2\ 

qs2 

/ 




a1/2 


(34) 


C(Ab) 

Ab 


+ (a tan a) 2 


Pb 


\Ph - Po 


4 

'c (Pb) " 

2 

■<^(Pco)' 

V. 

Poo/ 

_ Pb 

T 

_ P” 




2 


2 


2^ 

"a(a)‘ 

a 

+ 

~ 0 (q)‘ 
- *3 J 

+ 

■a(S)“ 
^ s _ 



Pb - Po 

s 

1/2 


(35) 


17 


a(L/D) = 


Fa 2 + Fn^ 


“i2 


(Fa cos a + Fjj sin a)^ 


a ( a) 2 + 


(Fa cos a + Fjj sin a) 2 


ct(Fn) 


% 

2 ^ 
a(FA)2l 

(Fa cos a + Fjj sin a) 2 


J 


1/2 


( 36 ) 


a(L/D) _ J 

(Fa^ + Fjq2 + 4FaFjj cos a sin a) a 

"a (a)' 

L/D ‘ ^ 

(Fjj cos ot - Fa sin a) (Fa cos a + Fjj sin a) _ 

a 


(% - 2Fa sin a cos cx)Fa 




cx>s a - sin a) (F;^ cos a + sin a) J 


L Fa J 


+ 


Fa% 

ta(FN)' 

2" 

_(F|j COS a - Fa sin a) (Fa cos a + Fjj sin ot) 

. Fn . 





1/2 


^(Cp) = Jqoo‘ 2 [a(P 2)]2 + q ^-2 a(pj 2 + 


(PZ - Poo) 




<^(Cp) 



( 37 ) 


( 38 ) 


( 39 ) 


Some of these formulations can be used immediately and in a general way. 
As an example, the National Transonic Facility (NTF) , a cryogenic tunnel using 
nitrogen is being built at this Center. Tentatively, the following parameter 
limits are being used for its design: 


1 8 


Specific heat of nitrogen Yn 2 atmosphere 1.4 to 1.45 

Specific heat of nitrogen at nine atmospheres 1.74 

Stagnation pressure p^ from M = 0.2 to 1.0, N/m^ 8.963 x 1 0^ 

Mach number 0.2 to 1.2 

R/A, m"'' 0 to 3.94 X 108 

q, N/m^ absolute 0 to 3.352 x 1 0^ 


The total temperature varies from 339 K to the temperature of liquid nitrogen 
for the value of Pt at that Mach number. 

Using the preceding information, programs were written in BASIC and run on 
the Hewlett-Packard 9830A calculator for 0 (M) , <7 (q) , and cr(R/i-) . The pro- 

grams are in appendix B and the plots of the standard deviations are shown in 
figures 1 to 4. These standard deviations are due to instrument errors only, 
and other independent sources contributing to the error of these variables must 
be included to obtain their total standard deviation. 


DESCRIPTICW OF INSTRUMENTATION 

Some of the basic instrumentation types whose parameters are used in this 
error analysis are described. The error values are real and were obtained 
from the manufacturer's error-band quotation and from data obtained from the 
instruments. 


The acoustic mancxneter is essentially a U-tube mercury mancmeter. (See 
fig. 5.) A sound signal is simultaneously sent through both legs of the manom- 
eter from which the differential mercury height, and thus the differential 
pressure, is obtained. This manometer is used over the following ranges, and 
the errors given by the manufacturers are shown as 3a values: 


1 . 034 X 1 q 5 N/m2 < 
4.826 X 1 q 4 N/m2 < 
2.758 X 105 N/m2 < 
1 .931 X 1 N/m2 < 


Pt < 2.758 X 1 q 5 N/m^; 
p < 2.758 X 105 N/m2; 
Pt < 4.826 X 1 0^ N/m^; 
p < 4.826 X 1 0^ N/m^; 


3a = 2.069 X 10l N/m2 
3a = 2.069 X lol N/m2 
3a = 4.137 X 10l N/m2 
3a = 4.137 X 10^ N/m2 


For higher pressures, a f used-quartz Bourdon tube gauge is used. (See fig. 6.) 
As the tube is deflected by a change in pressure, it is restored electromag- 
netically to a null position, thus minimizing errors. The current required to 
hold the tube in its null position is a measure of the pressure. Its 3a error 
is 0.006 percent of full scale plus 0.012 percent of the reading. In normal 
tunnel operations, for pressures greater than 4.826 x 10^ N/m^, only one fused- 
quartz Bourdon tube type of transducer is used whose full scale exceeds the 
estimated total tunnel pressure. In the NTF, however, two will be used with 
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full scales of 6.895 x 1 0^ N/m^ eind 1.034 x 1 0® N/m^. A sensor will activate 
a switch which will cut the lower pressure instrument off just before its full 
scale is reached. The data will be marked at the same time to identify the 
recording transducer. The full scales of the f used-quartz Bourdon tube pres- 
sure transducers start at 6.895 x 10^ N/m^, and increase in multiples of ten. 
There are also pressure transducers with full scales which start at 
1,034 X 10^ N/m^ and increase in multiples of ten. In fact, using three pres- 
sure capsules and staying within the pressure limits of these capsules, any 
full scale may be obtained by adjustment of the electronics. 


Acoustic mancmeters have at least three times the precision of the 
fused-quartz Bourdon tube type, but since their price is more than three 
times as high, they are used only when warranted by the precision require- 
ments. (Acoustic manometers have almost the accuracy and precision of a 
primary standard.) In tunnel work, acoustic manometers are used for measuring 
pressures below 4.826 x 1 0^ N/m^. For higher pressures, fused-quartz Bourdon 
tube pressure transducers are available with full scales of 5.171 x 10^ N/m^, 
6.895 X 10^ N/m^, and 1.034 x 1 0® N/m^. The last value is sufficient for 
tunnel work. However, these instruments can be obtained at very much higher 
full scales. By the same token, if the accuracy and precision of the fused- 
quartz transducer are adequate, they can also be obtained with full scales of 
6.895 X 1 q 4 N/m2, 1.034 x 10^ N/m^, 2.068 x 1 0^ N/m^, and 3.447 x 1 O^ N/m^. 
This is also essentially a static-type transducer. 


Stagnation temperature measurements are made with platinum resistance 
thermcsneters. (See fig. 7.) The operating range of these thermometers is 
from 8 to 533 K. In steady-state operation, a = 0.3 K. 


The six-ccsnponent strain-gauge-balance calibration fixture is shown in 
figure 8. This fixture, as well as manual loading, is used to calibrate the 
strain-gauge balances and obtain the coefficients for the interactions and 
second-order terms. 


RESULTS AND DISCUSSION 

This discussion is limited to some typical results rather than all the 
results that can be obtained and is used to illustrate what can be done. It is 
based on parametric data applicable to the NTF. In figure 1 , which is a plot 
of a(M) versus Mach number, 0(M) is dependent on the pressure parameter. By 
using basic relationships and dimensional analysis, and by following procedures 
shown in appendix A, a(M) can be written with dependency on other parameters. 
The sharp breaks in the curves of figures 1 and 2 occur because different types 
of pressure instruments, as well as the same type with different full scales, 
are used. Since the errors are a function of the full scale of the instruments 
used and/or their indicated value, a break occurs in the error curve when the 
pressure transducer is changed. Since these error calculations are based on a 
fixed total pressure, the breaks appear when the static pressure, which is a 
function of total pressure and Mach number, changes to a value requiring a 
change of transducer range with a corresponding change in error values. It is 
apparent that the greatest absolute errors in Mach number occur at the low 
Mach numbers. This is apparent from the relationships in equation (4) and 
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p = P|. (1 + and from the errors in p and p^. at different Mach 

numbers. It would be preferable not to try to interpolate values for pres- 
sures not plotted, because of the instrument error changes with pressure range. 
Instead, curves should be plotted for any additional total pressures of inter- 
est. The previous comments also apply to figure 2, which is a plot of a(q) 
versus Mach number for various total pressures. In figure 3, which is a plot 
of a(R/S.) versus Mach number at a temperature of 77.8 K, the greatest abso- 
lute errors are encountered at the highest Mach numbers and increase with 
higher total pressures. 

It is profitable to examine the standard-deviation equations in the body 
of the paper because they are algebraic. Discrete points of error buildup can 
thus be pinpointed to leurge multiplying coefficients. These points can be 
verified with actual numerical values. Time can then be more efficiently 
allotted for finding the means of eliminating or reducing these large sources 
of error. Examination of equations (13) to (18) shows that second-order terms 
as well as terms produced by interactions increase the errors in the strain- 
gauge-balance outputs. In addition, the errors made in measuring the interac- 
tion coefficients increase this error still more. Errors of this type result 
in nonlinear calibrations. There are other causes for producing nonlinear 
calibrations but these can usually be isolated. When nonlinearities attribu- 
table to the aforementioned causes are found, the production sources should be 
attenuated within the bounds of desired accuracy and precision and the limits 
imposed by the state of the art. 

Figure 4 shows that the standard deviation obtained for Mach numbers using 
only fused-quartz Bourdon tube pressure transducers of different full scales 
gives 2 to 3.5 times the error obtained in figure 1, where different types of 
pressure transducers were used. These larger errors still fall within accept- 
able limits for most research testing. 

Equations (1 3) and (24) to (28) give the relative errors for the strain- 
gauge-balance outputs. Although this may be the desired end information for 
some researchers, the actual loads and their errors are probably the more 
desired information. This information can be obtained from the calibration 
curves by finding the load corresponding to the balance output and applying 
the sigma values to obtain the extrema of the loads. The explicit, theoretical 
solution for the actual loads is omitted in this paper because it is believed 
that it more properly belongs in a paper dealing with the strain-gauge balances 
alone. These same equations show also that since the K coefficients appear 
everywhere as multipliers and are used in the iteration process, it is neces- 
sary that they be determined with a high degree of confidence in accuracy and 
precision. 


CONCLUDING REMARKS 

The work herein represents a beginning to a systenatic approach to error 
analysis for measurement instrumentation for the National Transonic Facility 
(NTF), which is located at the Langley Research Center. This work can be 
easily applied to any wind-tunnel measurements. It was long overdue in that a 
thorough search of the literature back to 1 934 revealed almost no such work. 


21 



Also, the work that was found constituted only a very small part of the effort 
being undertaken. Again, it is very simple to apply the error analysis of this 
paper to almost any instrument measurements. However, throughout this work 
randan errors are assumed; that is, there is just as much likelihood of the 
errors being negative as positive. They are variable in magnitude and dis- 
ordered in occurrence. It is also assumed that systematic errors are removed 
and good sampling practices employed. 

Pertinent wind-tunnel parameters, their standard deviations, and their 
theoretical derivation are given as are BASIC language computer programs and 
plots for obtaining the standard deviations of Mach number, dynamic pressure, 
and Reynolds number versus Mach number for the NTF. These standard deviations 
are for the instrument contribution only. It can be determined from the equa- 
tions and graphs which parameters are producing the largest errors or undesir- 
able errors and take the necessary steps for minimizing them or making them, 
insignificant. 

Finally, the author has noted that in many technical writings and investi- 
gations, the measurement errors are not reported, or are erroneously reported, 
partial errors are reported as total errors, or for whatever reason, the errors 
are glossed. The author hopes that this work, in some way, will encourage 
investigators to analyze their measurements, isolate and eliminate errors that 
have no role in statistical measurement, find the causes of and minimize their 
statistical errors when possible, and, finally, report their findings (their 
measured values) with their appropriate related errors. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 13, 1979 
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APPENDIX A 


DERIVATIONS OF STANDARD DEVIATIONS OF PERTINENT WIND-TUNNEL PARAMETERS 
Derivation of a(M) for All Mach Numbers 



(From ref. 4 , p. 4, eq. (44)) (Al) 
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From equation (7) 
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(All) 


Derivation of q Over Compressible and Incompressible Range 


q = 


YpM-' 


(From ref. 4, p. 4, eq. (31b)) 


(A12) 


In q = (In Y + In p + 2 In M) - In 2 


(A13) 


d(q) _ Sp 26 m 

q p M 


(A1 4a) 
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Derivations of Reynolds Number in U.S. Customary Units 



(A16) 


(Equation (AT 6) continued on next page) 
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+ 
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Tt 



lb 


J 


(A16) 


The only noncancelable dimensional item is p^/Tj.. if this item is in 
Ib/ft^/OR, then Reynolds number is in ft~^ . In the metric system, the deri- 
vation proceeds in the same manner. The 198.6 is a temperature in degrees 
Rankine and must be changed to degrees Kelvin to get a dimensionless ratio. 
This conversion yields 


R 

I 


1 . 2583 X 




+ 


198.6 

Tt 



ft 

Or lb 


1 

X (0.3048 m/ft) 

4.4482 N/lb 


(A17) 


Since reference 4 used the U.S. system and the constants contain degrees 
Rankine, when metric is used, an uncompensated T^. in the denominator is in 
degrees Kelvin and must be changed to degrees Rankine to correlate. Therefore, 


R 8.6221 X 104 

% 9 Or 

5 K~ 



1/2 


-4/7 


li 

(P 


-6/7 
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Tt 
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Derivation of Reynolds Number in SI Units 
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(Equation (Al 8) continued on next page) 
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(Equation (Al 9) continued on next paae^ 
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(Equation (A20) continued on next page) 
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Derivation of Ronaining Standard Deviation of Strain-Gauge Balance Outputs 
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+ <(Kg,3 + 2 FaK 6,9 + FuKgj4 


+ FyK 6,18 + MyKg^22 + MzKg^23 + %Kg^24)' 


_ 


^ra(FA)"^' 


Fa 


+ <(Kg^4 + 2MYKgjo + F^Kg^S + FYKg j 9 + FAKg^22 + MzKg ^23 




a (My) 


My 


+ ((^6,5 + 2MzKgji + FNKgjg + FYKg^20 


+ FAKg^23 + MyKe,25 + %K6,27) 


2rr 

Wj 


O (Mz) 


Mz 


+ /(I + 2MxKgj2 


FnK 6,17 FYKg^21 ■*■ FAKg^24 + %Fg^26 MzKg^27) 


Wi 


a(Mx) 


Mx 


(Equation (A25) continued on next page) 
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0(K6,1)‘ 

^ /FyK6,2\^ 

"cf(K6,2)‘ 

^ /FAK6,3f 

■^(K6,3)‘ 

. ^6rl 

^ \ Be / 

_K6,2 

\ Be / 

^6,3 


/MyK6,4\^ 

0(Ke,4)‘ 

^ /MzK6,s\^ 

0(K6,5)' 



V Be / 

. ^6,4 

\ Be / 

. K6,5 

^ \ Bg / 

.^6,7 . 


Vk6,8 \ 

^ Bg / 


jO(Kg,8)' 

^ /^FA%e,9\^ 

a (Kg ^9)' 


'‘^(K6,10)' 

1 

< 7 \ 

00 

1 

^ \ Bg / 

. ^6,9 

'’ \ Bg / 

. K6,10 . 


(^Z^^6,u\ 

'•^(Ke,!!)' 

fMx^K6,12^ 

<^(K6,12)' 

\ Bg / 

. *^6,n 

^ \ Bg / 

_K6,12 


rFjjPyKg J 3 


V 

0(K6,13)‘ 


"<^(Ke,14)‘ 

1 

_K6,13 . 

^ \ Bg / 

Kgj4 


2 


Pj^Y^ejs 

~6 



^P|^2Kg j 6^ 

‘‘^(K6,16)' 

. ’^6,15 

^ \ Bg / 

_Ke,ie 


2 


2 

/FnMxK6,1?\ 

*'^(Kgj7)1 

2 / \2 
^FyPaK6,18^ 

'<J(Ke,18)' 

\ Bg / 

.^6,17 J 

^ \ Bg / 

_K6,18 


^P^YK6,19f 

’o(K6,19)‘ 

2 2 
/FYMzKg^20^ 

'^(K6,20)' 

\ Bg / 

^K6,19 . 

^ \ Bg / 

K6,20 
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^^A%*^6,22^ 

CT(K6^22) 

\ ®6 / 

_K6,21 . 

\ Bg / 

- ^6,22 


/FaMzK6,23\ 

*^(^6,23) 

2a\ 

C^(K6,24)' 

\ ®6 / 

. ^6,23 

"^ \ Be / 

_Ke,24 


(n^zi^e,25\ 

“'^(K6,25)* 

^MyMxKg^2e'^ 

■^(K6,26)‘ 

\ Bg j 

. ^6,25 

"^ \ Bg / 

. Ke,26 


^M2MxKg,27 


.2 


vV2 


0(K6,27) 

_ ^6,27 J 


(A25) 
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TYPICAL BASIC PROGRAMS FOR OBTAINING SOME STANDARD DEVIATIONS FOR 


WIND-TUNNEL MEASUREMENTS 


1« PRINT "PROGRPlH FOR SIGMfl(H> VS.M AT VARIOUS STAGNATION PRESSURES USING" 

20 PRINT "DIFFERENT TYPES OF PRESSURE TRANSDUCERS. VALID FOR ALL HACH NUMBERS." 
SO PRINT "FOR USE ON THE HEWLITT-PACKARD MODEL 9830<;BASIC CALCULATOR)." 

40 PRINT 
50 PRINT 

60 PRINT " FILE 3" 

70 PRINT 

ftPH PPIHT "THIS PROGRAM WAS WRITTEN IN I.S.U. FEB. 28» 1979 BY E.RIND" 

90 PRINT "PROGRAM FOR SIGMA(M)VS.M VALID FOR ALL MACH NUMBERS." 

100 PRINT "Y=SIGMA<M>" 

110 PRINT "F=FULL SCALE ON PRESSURE GAUGE USED" 

120 PRINT "P1=ST AGNATION PRESSURE ;.N/Mt2" 

130 PRINT "P2=STATIC PRESSURE) N/Mt2" 

140 PRINT "S1=SIGMA<P1)5N/Mt2" 

1 50 PR I NT " S2=S I GMA ( P2 ) ; N/Mt2 " 

160 PRINT "G=GAMMR" 

170 SCRLE 1.2J0J3.5E-04 

1 8 0 R = 4 . oil 6 E + 0 5 

190 B=6.895E+05 

2 0 0 C = 1 . 0 3 4 E + 0 6 

210 n=2.758E+05 

220 H=1.931E+015 

230 J=1 . 034E+05 

240 K=4.S26E+04 

2 5 0 N = 2 . 0 6 SE + 0 5 

260 G=1.4 

270 Gl=l/G 

280 RERD PI 

290 PRINT 

300 PRINT ‘■pl==’’Pr’N/r-1t2;STRGNRTI0N PRESSURE" 

310 PRINT 

320 FOR N=0.2 TO 1.2 STEP 0.1 
330 P2==Pl*a+M*M/5)t(-7/2> 

340 IF P1>R RND P2>R RND PI <= B RND P2 <= B THEN 380 

350 IF P 1 >B OR P2>e THEN 360 

360 F=C 

370 GOTO 390 

o y 0 F = B 

390 S2=(6E”05*F+P2*1 . 2E-04>/3 
400 S1 = (6E”05*F + P1*1 . 2E“04::'/3 
410 IF P1>R THEN 500 

420 IF PI <= R RND P2 <> R RND PI >= D RND P2 >= H THEN 480 
430 IF PKD HMD P2 <= D RND PI >= J RND P2 >= K THEN 450 
440 IF PKJ OR P2<K THEN 450 
450 S 2=6. 895 
460 S 1 = 6 . 8 9 5 

470 IF S2 OR SI =6. 895 THEN 500 
480 S2=13.79 
490 SI =13. 79 

500 Y=(Pl/P2)t(*-Gl )*(1..-'<M*G*P2> ) *SQR < S 1 12+ (P 1 /P2 ) t2*S2t2 > 

510 PRINT ’■r'i="M; "y="y; "pi="pi; •■p2="P2; "gi = "gi; ■•g="Gj "si = "si5 ■■s2="S2 

520 GOTO -540 

530 PLOT MsY 

540 NEXT M 

550 PEN 

560 DRTR 1 . 035E+05 j 2. 068E+05n 2. 758E+051. 3 . 447E + 05^ 4. 137E+05 ^ 4 . 826E+05 ^ 6 . 205E+O 
570 DRTR 6 . 895E+05 ? 8 . 274E+05 ? 8 . 963E + 05 ? 1 . 034E+06 
580 GOTO 280 
590 END 

600 PRINT "THIS PROGRRM NRS RUN BY 


ON 
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10 PRINT "PROGRAM FOR SIGMA CM) VS. M FOR DIFFERENT STAGNATION PRESSURES" 

20 PRINT "USING VARIOUS FULL SCALE RUSKA PRESSURE TRANSDUCERS , VAL I D FOR ALL MACH 
30 PRINT "NOS. FOR USE OH THE HEWLITT-PACKARD MODEL 9830 (BASIC:) CALCULATOR." 

40 PRINT 
50 PRINT 

60 PRINT " FILE 7" 

70 PRINT 

80 PRINT "THIS PROGRAM WAS WRITTEN BY E.RIND MAY 17» 1979, " 

90 PRINT "PROGRAM FOR SIGMA(M:)VS. M USING VARIOUS FULL SCALE RUSKA PRESSURE" 

100 PRINT "TRANSDUCERS. VALID FOR ALL VALUES OF MACH NUMBER." 

110 PRINT "ALL UNITS ARE IN I.S.U." 

120 PRINT "Y=SIGMA(M:)" 

130 PRINT "Pl=STAGNATION PRESSURE-N/Mt2" 

140 PRINT "P2=STATIC PRESSURE-N/Mt2 " 

150 PRINT "Sl=SIGMA(Pi;)" 

160 PRINT "S2=SIGMA(P2:)" 

170 PRINT "M=MACH NUMBER" 

180 PRINT "G=GRMMA=CSUBP/CSUBV" 

1 90 SCALE 0 » 1 . 2 ! 0 j 0 . 00 1 33 

200 G=1.4 

210 Gl=l/G 

220 READ P1»F,F1 

230 PRINT 

240 PRINT 

250 PR I NT " P 1 = " P 1 " ST AGNAT I ON PRESS . -H/MT2 " !• " F= " F j " F 1 = " F 1 


260 PRINT 

270 FOR M=0.2 TO 1.2 STEP 0,2 
280 P2=P 1 * ( 1 +M*M/5 ;■ t < -7/2 

290 DATA 8.963E+05! 1.034E+06? 1 . 034E+06 j 8 . 963E+05 j 1 , 034E+06j 6 . 895E+05 
3100 DATA 8. 963IE+05? 1 , 034E+06? 3. 447E‘+'05 

310 IF Pl=8.963E+05 AND F=1.034E+06 AND Fl=1.034E+06 THEN 590 

320 IF Pl=8.963E+05 AND F=1.034E+06 AND Fl=6.895E+05 THEN 610 

330 IF Pl=8.963E+05 AND F=1.034E+06 AND Fl=3.447E+05 THEN 640 

340 REM THE TOTAL PRESSURE IN THE TUNNEL FOR THE FOLLOWING IS 5.171E5 N/Mt2 ABS. 

350 DATA 5 . 171E+05 > 1 . 034E+06 j 6. 395E+05j 5 . 171E+05? 6 . 895E+05 ? 6 . 895E+05 

360 DATA 5. 171E+05» 6. 895E+05> 3. 447E+05j 5. 171E+05> 6. 895E+055 2. 068E+05 

370 IF Pl=5,171E+05 AND F=1.034E+06 AND Fl=6.895E+05 THEN 610 

380 IF Pl=5.171E+05 AND F=6.895E+05 AND Fl=6.895E+05 THEN 670 

390 IF Pl=5.171E+05 AND F=6.895E+05 AND Fl=3.447E+05 THEN 700 

400 IF Pl=5.171E+05 AND F=6.895E+05 AND Fl=2.068E+05 THEN 730 

410 REM THE TOTAL PRESSURE IN THE TUNNEL FOR THE FOLLOWING IS 2.068E5 ti Mt2 ABS. 
420 DATA 2.068E+05 j 1 . 034E+06 j 3 . 447E+05) 2 . 068E+05> 6 . 895E+05 j 2 . 068E+05 
4310 DATA 2 . 068E+05 !< 6 , 895E+05 1 1 . 034E+05 

440 DATA 2.068E+05! 3.447E+05) 2. 068E+05j 2, 068E+05J 3. 447E+05j 1 , 034E+05 

450 DATA 2. 068E+05- 2. 068E+05? 2. 068E+05> 2. 068E+05! 2, 068E+05? 1 . 034E+05 

460 IF Pl=2.068E+05 AND F=1.034E+06 AND Fl=3.447E+05 THEN 640 

470 IF Pl=2.068E+05 AND F=6.895E+05 AND Fl=2.068E+05 THEN 730 

480 IF Pl=2.068E+05 AND F=6.895E+05 AND Fl=1.034E+05 THEN 750 

490 IF Pl=2.068E+05 AND F=3.447E+05 AND Fl=2.068E+05 THEN 790 

500 IF Pl=2.068E+05 AND F=3.447E+05 AND Fl=1.034E+05 THEN 820 

510 IF Pl=2.068E+05 AND F=2.068E+05 AND Fl=2.068E+05 THEN 850 

520 IF Pl=2.068E+05 AND F=2.068E+05 AND Fl=1.034E+05 THEN 870 

530 REM THE TOTAL PRESSURE IN THE TUNNEL FOR THE FOLLOWING IS 1.034E5 N- Mt2 ABS. 
540 DATA 1 .034E+05! 2.068E+05»2.068E+05) 1 . 034E+05 j 2. 068E+05> 1 . 034E+05 
550 DATA 1.034E+05J 1.034E+05ii 1.034E+05 

560 IF Pl=1.034E+05 AND F=2.068E+05 AND Fl=2.068E+05 THEN 850 

570 IF Pl=1.034E+05 AND F=2.068E+05 AND Fl=1.034E+05 THEN 870 

580 IF Pl=l,034E+05 AND F=1.034E+05 AND Fl=1.034E+05 THEN 900 

590 81=82=137.9/2 
600 GOTO 910 
610 S 1=137. 9/2 


620 32=96.53/2 
630 GOTO 910 
640 Sl=137.9/2 
650 S2=55.16/2 
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FILE 7 CONTINUED 

..I,. 1 ,1 1 { n 'MM 


.‘H 

[ 

' h . 

r.; 

7 i f i 1 

i l•■1 

•-IU 

TOP 



7 

b : 

? 1 0 

. _ L 1 1. 1 


720 

M 1 [n ■ 

qn 

?:"i0 

.J 


74U 

- ”:■! « 

'l3 

7‘.n 

i 1,1 

■i ! 0 

7bU 

.i „ 

! 1 . 

7.4j 

■ . ' - . '0., 

be;! 

780 

.Mill 

1 0 

790 


, 1 6 

800 

8 2 "" 3 4 . 

, 48 

810 

GOTO 

:j 1 0 


840 GOTO 910 
860 GOTO 910 

890 GOTO 910 
9 0 0 S 1 = S 2 = 2 0 . 6 8 .■■■' 2 

9 1 0 V = ( P 1 . ■ p 2 ) -t' ( G 1 ) < 1 ( M * G * P 2 > > * S Q R < S 1 -t- 2 + ( P 1 . - P 2 > 'T' ? * 8 ? 1 2 > 

920 PR I NT " N= " M j " Y= “ Y i " P 1 = " P 1 ; " F= " F ; " P2= " P2 j " F 1 = " F 1 1 " S 1 = ’' 5 1 - " 82= " 8? ; 

930 PRINT "G1=“G1 ; "G="G 

940 PLOT MjY 

950 NEXT IT 

960 PEN 

970 GOTO 220 

980 PRINT 

990 END 

1000 PR INI THIS PROGRfiM NflS RUN BY ON 
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Ifi PRINT " THIS IS THE PRQGRRM FOR SIGHR (Q>VS.t1 FOR DIFFERENT STflGHftTlQU PRESS- 
20 PRINT "ORES AND DIFFERENT TYPhb UP PKhsyUKh IRHNSDUCERS. VfiLID ALL MfiCW" 

30 PRINT ' tJIJHUtRT-:'. IN i.y.u. antiKll SVyiEH). HJk USE UN I Ht HENLE I T-PHCKflRIi" 

4 0 p R I N T ' HuIiEL 0330 (.tiHSU; .:' uhLuulh I uk . ' 

50 PRINT 
60 PRINT 

70 PRINT " FILE 8" 

80 PRINT 

90 PRINT "PROGRAM FOR SIGMfl(Q) VS. M IN I.S.U.BY E.RIND flUG15jl979" 

100 PRINT “USING FORMULAS APPLICABLE OVER ALL MACH NUMBERS" 

110 PRINT 

120 PRINT "Y==SIGMA CQ)" 

130 PRINT "F=FUL.L SCALE ON PRESSURE GAUGE USED" 

140 PRINT "Pl=STRGNATIOH PRESSURE, N/Mt2" 

150 PRINT. "P2=STATIC PRESSURE, N/Mt2" 

1 6 0 P R I N T " S 1 = S I G M A ( p 1 > , N / M 1 2 " 

1 7 0 P R I N-T " S 2 == S I G M A ( P 2 ) , H / M T 2 " 

180 PRINT "S3=SIGMA(;M) " 

190 PRINT "Q=THE DYNAMIC PRESSURE, N.'Mt2" 

200 SCRLE 
210 . 034E + 05 

220 B=2.758E+05 
230 C-3.447E+05 
240 ;[i 4 . 0 2 6 El +' 0 5 
250 S95E + 05 

2 b 0 J 8 . 9 b 3 E •+ 0 5 
270 K--1.034E + 06 
280 N-1 .931E+-05 
290 lj:“4 , 826E + 04 
300 

310 i::7E=^^^i - -G 
320 REEili PI 
3.-:0 PPIN! 

:M0 PE-iN I ''pi-''pr'H/r-it2; strgnhtion pr:essure" 

3SM ppm! ■3:;Y^Pl?P2.Sl,S2mNIi S3 RE;:E in 
>;. 0 pE’[Mf 

■::,^0 Mip N-I1„2 TO 1.2 STEP 0.5 
•::PM P2:-Pl l fN'^M/5::'t(-~-7/2> 

.TH0 2‘P2*Mt2 

400 IF Pi>D RND P2>D RND PI <= H RNH P2 <> H THEN 440 

410 IF P1>H RNIi P2>H THEN 420 

420 F-K 

430 GOTO 450 

440 F-H 

4 5 0 2 ™ < 6 E 0 5 F + P 2 * 1 . 2 E - 0 4 ) /' 3 

4 b 0 1 - C E 0 5 F + P U l . 2 E - 0 4 ) / 3 

470 IF P1>H THEN 570 

480 IF PI <= II RNIi P2 <= D RND PI >== B RND P2 >== N THEN 550 

490 IF Pl:B RND P2<B AND PI >= R RND P2 >= U THEN 510 

500 IF PKR OR P2<IJ THEN 510 
510 S2^^^b.895 
520 SI -b.R3b 

530 IF '..■•b.bMS then 570 
540 U M then 570 

550 S2-.1. P :'3 3 
5b0 M 1 r::-' 

5 70 S ! • • P 1 P 2 - 1 1 - G 1 ) ± ( i / ( M ± G * P 2 > > 3 Q R ( S 1 1 2 •+• < p 1 .••• p 2 > -t' 2 S 2 2 > 

^-1 y 0 r ' < I 4' ' ' ' ‘ I j ,J I I’l ’t' v ;i T' 2 •+• b E' 2 M o 3 T 2 

530 PLU I N:i Y 

000 PP'INl "I'P ‘N'N "Y^^"Y; "P1-"P1; '•P2-“P2; “S1 = ‘’S1^ “S2==:=“S2 i “S3===3'S3 
bl0 PEN 
F20 MEYT M 

i-Y : 0 I) R 1" 1 . 0 ::: 4 E ■+■ O ^5 1. 2: . 7 5 8 E + 0 5 ? 3 . 4 7 7 E •+• 0 5 ? 4 . 8 2 b E •+• 0 5 ? b . 8 9 5 E •+■ 0 5 ? 8 .. 9 6 3 [■ •+■ 0 . 1 .. 0 3 4 E + 0 6 
i.40 GOTO 320 
!-.50 END 

h.,G0 pRIlPr 'THIS PROGRRN MRS RUN BY 


ON 
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THIS IS THE PROGRHil FOR 
5 T H G H fl ! 1 0 H PRESSURES U S 
vHLID HJK ALL HflCH HU rib 
aiHSlC.) CHLCULHTOR. " 


IiIFFEREHT 


II I 

OF PITS' 
HI WLE I r 


:o PRINT 

'0 PRINT '■ FILE 5" 

;0 PRINT • 

,0 PR I HI ''III I', F'Pi ii ,|<FiM HHS NR I I TEN BY E.RINIi I'lHRCH 

00 PRINT "l'l■■llU■■llll FOR SIU'IH'F' L ) VS „ MflCH NO . IN I.SJJ. FOR fill VRIHFS uF Mfir.l 
10 PRINT 'Tl S nil, NOTION F'F'Fl:-.Sl IRE < N/Mt 2 " 

2 0 P R I N I “ I' 2 -• S T H T I C P R E S S U R t - N / M 1 2 " 

30 PRINl "■■.J--SINI'1H<P1>,N. Mfj" 

4 0 P R I N I '■ 2 S 1 0. M 0 ( P 2 ) j N / M 1 2 " 

50 PRINl "M-. [iVNOMIC PRESSURE? N/Mt2" 

60 PR I N I " I r'UBSOLUTE TOTAL TEMPEROTLIRE , ( DEG . KELV I N ) " 

7'M Fi:™ 1 . !■! :-l IF i 
::n i: . ' . .‘S! :L -i U‘j 
!''H I •i-ri-'L + uS 

kl k ! i ! ‘ I „ I" ! i;"i "I" k 1 l“i 

II U J 

'.-■I-J M -IL ^:*b3E;' + yS 

'^V|Hu 5 


■•jn I. 

"™4. 

l;; i E ■+■ k1 4 







kU 1” 

1<IHT 

■■RLL TEHF- 

ERFITURE 

3 MUST 

BE.. 

IH BEG 

. KEl 

,’M F 


" ' {-'6 I b!''lH 







::n E 

:EH Tl 

78 HE 

G. F<ELV 

I FL== 

140 

HEG 

.. RFiNK 

I NE 

■lu F 

TIM 1. 

I J ! . I HE 

G. F(ELV 

I 

200 

HEG 

. Rh\-W- 

INE 

;nn F 

:e:m [ 

HE 

G. F<ELV 

I 

300 

IJ E I..1 

.. RRNh: 

I NE 

: 1 H F 

ti;m n 

.■‘.•2.2 HE 

G. F;:elv 

I N=:= 

400 

HEG 

.. RRNK 

I NE 

>20 F- 

■■F:M I‘| 

I .•7 7.:': HE 

G. FaiLV 


000 

HEG 

.. RRNK 

I NE 

0 i- 

■TiiM h. 

. HE 

G. FCELV 

I 

600 

HtIG 

. RRNK 

I NE 

;:40 1' 

■FM 1. 

■ • 1 BE 

G. KELV 


610 

BEG 

RRNK 

:INE 

;:50 F 

IN 1 

"F - 1 t.ll.L S 

CRLE UN F-'R 

ESS! 

IRE 

I..J Ii 1 J 1.1 E 

USED 


6 0 R EM ' R' I ' 1,2 5 8 3 E 0 6 * P 1 * M / T 1 2 * < 1 + M * M / 5 > t < - 1 , 5 ) + ( ( 1 + M * M / 5 ) t ■: - 1 * 1 + 1 '4 fi , 0 1 

7 0 R E M ( R / L ;:■ = 1 . 2 5 8 3 E 0 6 * S Q R ( ( ( P 1 ' P 2 1 ( - 4 / 7 ) - i P 1 P 2 ) t - 6 / 7 ) ) * 5 I 

8 0 R E M * ( P 1 1 ( 5 7 ) P 2 1 ( 2 / 7 > . ■ ■ ■ T + 1 9 8 , 6 * P 1 / T 1 2 

3 9 0 P R I N T " S 1 = S I G M fl ( p 1 > 5 S 2 = 8 1 G M Fl ( P 2 ) ? S 3 = S I G M fl ( T ;■ TI = T E M P ? D E G . K E L V I N ? M == l'1 fl i ; H N 0 

400 PRINT "G=GFlMr'1fl" 

410 REM S3 =0.5 riEG./3(DEG.C)=0.9/3 DEG. F=0.3 DEG, RFlNKINE 
420 G=1.4 
430 Gl=l/G 

440 SCALE 0 j 1 . 2 ) 0 j 4 . 9E+06 
450 G2=G1--1 

460 PRINT "W=REYNOLD’S NUMBER.-'L" 

470 FOR J=1 TO 7 
480 READ T 
490 FOR 1=1 TO 7 
500 READ PI 
510 PRINT 

5 2 0 P R I N T " p 1 = " P 1 '■ N / M 1 2 ? S T fl G N fl T 1 0 N P R E S S U R E " ? " T = " T " T E M P ? D E G . K E L V I N " 

530 PRINT 

540 FOR M = 0.2 TO 1.3 STFP O F, 
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FILE 5 CONTINUED 

' 1 I ’ j.- 1 ( ;i .f|.| + |.| ' , I t ( > 

* ii .u I !• ' I "I M '•!■! ' I' ' 7 •••■■2 ••••F'u' - F 1 P-! 

IF IM D Rlili 1‘2 D FlHB PI H RHD P2 <= H THEH 616 

i| IM HRN.IJ I H THEN 
M-iu !• ■•■! 
t.;! 6 ! •, I ,! I it |/i 

F-H 

I • u ■ R I'M M R M M ' . ' 1 , 2 E ” 6 4 > •••■■ 3 

i ' ' i"T' + F +1 i •- 1 , 2E”04 ) •••■'3 

1.4R IF PI I) RNIi P2 n RNIi PI B RND P2 N THEN 766 

i! !•] B MNI) P2 < == B RNI) P1>R RND P2>U THEN 676 

i.^.u !i FM R Uh: P2<IJ THEN 676 

!•••. ,• ' M !" i ] i-i „ I 'l ' "I !.; I 

6’. ‘.M 3 i 1 . .. 

hi.RR :M ! UP F. 2 6. 5 THEN 728 

i’Ou 1 :„739 
.'|M .I:-! 789 

.•2u ■ . M.. lb,' 

F . n i I U ! M . I ' n ' u 4 J 9 3 4 E + 0 5 ? 2 „ F 5 8 E. + 9 9 :* • .. 4 4 F E + 0 5 ? 4 , :j . • F i 0 5 ? 6 „ y 9 5 E + 0 5 k 8 == 3 6 • ii: FI F:i 5 
7-1 M iUi III 1. - u- :4 Eh-06 

.’RU iiU ! H 1 1 1 « 1 - 1 . 034E+0R!t 2. 7R8E+0 Rh 3.. 447E + 65 ? 4 . x.b.F ^ 65^ 6 . S9RE+05 M4 . 963E+05 
[1)1 II- 1 I. .. ki ^FUlMTJb 

i ' u J j l -l I H I F. i "; „ 1 I. „ i;;i :l! 4- E + 9 5 ? 2 . i-'' 5 8 E -i- 9 l::;i !< : Is .. 4 4 7 E -i- 9 5 ? 4 » : : i- • F i 9 R h f , . y 9 R E + 9 5 :> 8 .= 9 6 :3 E + 9 5 

.'uu iHi!l-l lMi.:i34E + 06 

F u n it)! "1 \ 1 2 3 „ 3 , 1 .. 0 3 4' E + 0 R ? 2 .. 7 5 8 E + 9 5 :< 3 „ -4 4-7 E •+• 05 h 4 « S 2 G E "i- 9 5 !< 6 „ 8 9 5 E -i- 9 5 n 8 .. 9 6 .;:i E. + 9 5 
uiiu i!ii 1 11 1. .. i::i34E K;i6 

• i- u i t l-l i l"l I: ■“ i"' . ' 1 « R 4- E +• 9 5 1 2 . I’ i;-i E’ -i FJ .. 44 1’’ E -i- 9 5 ? 4 « 31.' 6 E - i- k;i 5 •< G .. 8 ‘‘.h R E + 9 5 ? i;..i . G ' i F- -i u 5 
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Figure l.“ cr (M) versus Mach number for various stagnation pressures using different 

types of pressure transducers. 
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Figure 2.- a(q) versus Mach number for various stagnation pressures using different 

types of pressure transducers. 
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Figure 3.- a(R/S,) versus Mach number at T = 77.78 K (140° R) for various stagnation 
pressures using different types of pressure transducers- 
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Figure 4.- Concluded 







Figure 5.— Acoustic U— tube mercury pressure manometer 
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Platinum resistance thermometer with sheath 






component strain-gauge-balance fixture being loaded. 
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